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Дискретные аналоги интегральных 
уравнений Вольтерра П рода 


Предлагается класс математических моделей процессов в линейных системах с фиксированными 
состояниями. Такие математические модели являются дискретными аналогами интегральных уравнений 
Вольтерра П рода. Рассмотрены методы аналитического решения предложенных уравнений на основе 
отыскания решений соответствующих уравнений, связывающих ядро и резольвенту. Рассмотрена и 
доказана теорема, устанавливающая аналитический вид резольвенты по заданному ядру для ряда 
важных частных случаев, а именно при сепарабельном виде ядра. Рассмотрены эквивалентные 
преобразования предложенных математических моделей к линейным разностным уравнениям. 


Введение 


Проблема исследования управляемого изменения состояния сложных динами- 
ческих систем в настоящее время решается путем построения их математических 
моделей (ММ), как правило, в виде моделей пространства состояний, как отвеча- 
ющих в наибольшей степени методам и средствам современной теории управления. 
Однако далеко не все процессы в реальных объектах управления могут быть описаны 
указанными математическими моделями пространства состояния (ММПС), что требует 
отыскания новых форм математического описания их движения. 

Важная научно-прикладная задача состоит в отыскании таких форм ММ, которые 
при сохранении адекватности реальным процессам позволили бы осуществить их 
численную реализацию. 

Известные преимущества интегральных моделей, в частности интегральных урав- 
нений Вольтерра П рода, предопределяет необходимость и практическую значимость 
рассмотрения методов отыскания различных форм их аналитических решений, на 
основе решения соответствующих уравнений, связывающих ядро и резольвенту. 


1 Постановка проблемы и цель исследования 


Теоретические основы решения интегральных уравнений рассмотрены в ряде 
фундаментальных работ [1-3]. Методы и алгоритмы их численного решения пред- 
ложены в справочной литературе [3], [4]. Особенности интегро-дифференциальных 
уравнений и программные средства их решения изложены в [5], [6]. Для отдельных 
частных случаев интегральных уравнений в [1-3] предлагаются также некоторые 
аналитические решения. В [7], [8] предложены такие решения для ряда важных 
прикладных задач. 

В то же время необходимость численной реализации ММ непосредственно в 
составе систем управления требует систематического рассмотрения вопросов отыс- 
кания решений интегральных уравнений и их дискретных аналогов, применяемых 
в качестве математических моделей исследуемых объектов. В первую очередь это 
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касается интегральных уравнений Вольтерра П рода, имеющих широкую область 
применения в прикладных задачах и для которых разработаны эффективные методы 
численного решения [3]. 

Целью настоящей работы является разработка методов аналитического решения 
дискретных аналогов интегральных уравнений Вольтерра П рода на основе отыскания 
решений соответствующих уравнений, связывающих ядро и резольвенту. 


2 Основные результаты исследования 


Интегральное уравнение Вольтерра П рода 
у(х) = Х(а)+ [К (5.5) 55) ( 
ное СО и оно 
ое С, © =сотзь 


К (х,5) — ядро, непрерывное в пределах треугольника х> а,5 <В,х=5.Ь>а, 
х,5 — независимые переменные, а = сопз6, имеет следующее решение в анали- 
тическом виде 


у(х) = бо [быль (2) 


гдег (х, 5) — резольвента, удовлетворяющая уравнению 
х 
г (х,5) = К (х,5) + [ (хр) т (1,5) 4, (3) 
а 

Известные [1-3] условия существования и единственности решения уравне- 
ния (1) заключаются в непрерывности ядра К(х,5) на сторонах и внутри треугольника 
х>а,5 <Б,х=5.Ь > а, а также в условиях ограниченности функции /(). 

Решение уравнения для резольвенты (3) имеет важное значение как для теории 
интегральных уравнений, так и для прикладных задач, поскольку отыскание резоль- 
венты по известному ядру из (3) предоставляет возможность решить тем самым целый 
класс интегральных уравнений Вольтерра П рода с различными правыми частями в 
виде функций Дх). 

Рассмотрим дискретный аналог уравнения (1) 


Хп 
у(х,) = Л(х,) +. К (х,,5, (5, }, (4) 
м 

где п] Е М, 1 < пи. 

Уравнение (4) может рассматриваться, с одной стороны, как приближенная квад- 
ратурная реализация уравнения (1). С другой стороны, (4) может использоваться как 
математическая модель процессов в системах с дискретными состояниями, например, 
в цифровых системах управления, моделирования потока отказов, отражения сигнала 
от дискретных рассеивателей и др. В первом случае нахождение аналитического 
решения (4) позволяет избежать известных проблем численного решения уравнений 
с переменным верхним пределом. Во втором случае нахождение решений уравнений 
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вида (4) имеет важное самостоятельное значение, поскольку дает возможность по- 
строить алгоритмы нерекурсивного типа, соответствующие требованиям реального 
времени. В последующем (4) будет определяться как % -уравнение. 

Будем полагать далее, что справедливы следующие условия 


У /(х) <С,, С, =сопЯ (5) 
"(х„, 5," (х, 15 55 С., С; =соп& (6) 
пы С. = сопя (7) 


В дальнейшем переменные-аргументы опускаются и используются следующие 
обозначения 
и ‚=, ый. 


У(х,)=у,, Л(х,) = Л, 


Решение уравнения (1) будем отыскивать по аналогии с (2) в форме 


РИ (8) 
= 
Из (4) и (8) следует уравнение для резольвенты 
= №" (9) 


=] 
Вывод уравнения (9) может быть получен следующим образом 


Хил) ХК, к ‚= 


И 
м 
= 
— 
Е 
ие 
= 
— 
м 
о 
ГИ: 
— 
п 
= 
т 
ее 
ах 
592 
ы 
— 


х,.5)) (5, м ыы 


5,=Х 5, = У; =5) 


| 
м 
= 
— 


С учетом (5) из последнего соотношения следует (9). 

Для исследования свойств и нахождения решений >» -уравнения типа Вольтерра 
рассмотрим первоначально следующую Лемму. 

Лемма. Для однородного »-уравнения Вольтерра П рода решением является 
0-последовательность. 

Действительно, для уравнения 


— >. КУ, 
р 
решение есть 
у, =0, УпеЕМ, 


что непосредственно следует из (8) при Г, =0, УпЕМ. 


Решения уравнения для резольвенты (9) устанавливает следующая теорема. 
Теорема. Если ядро является разделяющимся (сепарабельным) и имеет вид 


и. (11) 
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где @, =со(а, 3 т Ее аб ный] и справедливы условия 


"> “ит 


(5 —7), то решением уравнения для резольвенты является решетчатая функция 


а Е 0 (12) 
где матрица Е, определяется из соотношения 
п ое 

Е, = [(Е-Аа’). (13) 


=) 
Необходимость. Прямой подстановкой (12) в (9) получаем 
Е, Е +В.аЕ,, (14) 


Откуда непосредственно следует (12) и (13). 
Достаточность. Пусть (9) имеет два решения 


Ш 
(0 — (0 
7; ЕК р. в 


т й 
=] 


тогда 


РО —-"® = Ан, = У Е,Аи,. (15) 


п] п] п] 
Уравнение (15) является однородным уравнением и согласно Лемме имеет нуле- 
вое решение по координате х,. 
Рассмотрим решение по координате 5,. Пусть п =1, тогда Ди, = К,Аг,, отсюда 
Ак, =0. 


Далее получаем 


Ах, = АГ, АЕ 
ДУ, 1 7 Ду, 7 т А „Ат, 1 7 


В предположении, что Ду, =0, следует Ду, ,, =0. По индукции устанавливаем 


+1, 7 
Ду,, =0. Следовательно, решение (9) единственно и определяется (12) и (13). 

Следствие. Если К, и 7, — ядро и резольвента уравнения (9), то А„, =фФ_К,„ и 
7; =ф,К„, также являются решением (9) при ф, #0, УпЕМ. 


Рассмотрим > -уравнение Вольтерра П рода с разделяющимся ядром: 


у, =/,+а, >». Ву, (16) 
Обозначая 
=. ВУ, (17) 
= 


получаем следующее разностное уравнение, эквивалентное (16) при нулевых началь- 
ных условиях 


У, = В, (ал, + Л, )+ у, (18) 
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Из (16) и (18) следует дискретная математическая модель 


У, = И ВЕЕТ Ри ® 
1-@, д, 1-а, р, (19) 
У» == ау, ть ие 
с решением 
=Л+а,> ЕВ, (20) 
= 
где 
й —1 
Е, = ](1-ар,) . (21) 
=] 


Рассмотрим более сложное, по сравнению с (16), У -уравнение Вольтерра П рода 
с разделяющимся ядром 


у,=.+1 У Ву, (22) 
= 
Обозначим 
7, =>. Ву, 
= 


Отсюда следует дискретная математическая модель пространства состояний В виде 
- —1 - Г - 
ор -т\ | - т 
у, =(Е-В.а!) у, +(Е-Р,а!) В 
= -т- 
У» — Ч м, т Л * 
фундаментальная матрица которой имеет вид 
- —1 Г 
г. ЕР т. 
Ф„=(Е-Р,а, > Е, =] Ф.. 
1=] 
Решение системы (24) определяется соотношением 
п 
=: РАЙ 72. 
ЕЙ ИВР: 
1: 


В силу эквивалентности (22) и (24), это решение одновременно является реше- 
нием > -уравнение Вольтерра П рода с разделяющимся ядром. 

Возможности предлагаемого подхода иллюстрируют следующие примеры. 

Пример 1. Пусть ядро в уравнении (16) задано в виде 


(24) 


р 
А =——а "а'. 
у 
а 
Используя результат доказанной теоремы, получаем резольвенту 
„=а-1. 


При задании «правой части» следующей решетчатой функцией 
—_ —п 
ДЕЁ-Г 
получаем решение в виде: у, =1, УКе М. 


Настоящий пример является дискретным аналогом уравнения, рассмотренного 
в [3] стем же решением, полученным методом квадратур. 
Пример 2. Пусть ядро в уравнении (16) задано в виде 


к, =-(а-а "а7 = (1-а)а 2. 


пт] 
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На основании установленных аналитических результатов получаем резольвенту 
Е. =а" 1 =а®).а” 


2п 2 —2(п-7 
г, =(1-а)а "а" /а=—— И 
у 
а 
При задании «правой части» следующей решетчатой функцией 
Г. =К-а”" 


получаем следующее решение: = а" . 
У» 


Пример 3. В условиях примера 2 
К. = (1-а)а"а’ ь 


п} 
г, =(1-а)а а” Га, 


зададим другой вид «правой части»: {, =1. В этом случае решение имеет вид: 


Е 


Последний пример иллюстрирует важное достоинство резольвентного решения: 
могут быть получены решения классов У-интегральных уравнений с различными «пра- 
выми частями». На рис. 1 и рис. 2 представлены ядро и резольвента для а = 2. 


_2-П 291 


-2000 - 


-4000 - 


-6000 


-8000 


5 5 


п 9: 


Рисунок 1 — Разностное ядро для примера 3 


22 1+2 9:1 


5 о 


п о 


Рисунок 2 — Резольвента для примера 3 
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Если К, =К,_, — разностное ядро, то в силу линейности Х-интегрального уравне- 
ния т, =/,_, — резольвента также является разностной. Отсюда уравнение для резоль- 
венты имеет вид 

Ш 
И р а. (25) 
=] 

Полагая, что г, и К, удовлетворяют условиям существования их 7-преобразо- 
ваний, где 5=и-), из (25) следует операторное уравнение: 

В(2)=К(2)+К(2)В(2), (26) 
где К(2),К(2) — 2-изображения ядра и резольвенты соответственно. Отсюда 


следует аналитическое решение уравнения для резольвенты в виде 


в=7'{В(2)}, (27) 
где 
к (=) Е 8) 


Решение (27) и (28) обобщают операционный метод решения интегрального 
уравнения Вольтерра П рода с разностным ядром на соответствующие У-уравнения. 
Пример 4. Пусть 


КИЕК,,=а, ТЕГ, =и(1-а)"”. (29) 


21 п 
Определим изображение ядра 
” 
К(2)=а—^_. 
2-1 
Согласно (28) получим изображение резольвенты 


а2-1 а 2 
Ю — в : 
(:) (1-а)2-1 1-а Е 
1-а 


Используя обратное преобразование, получим окончательный результат 


о (Е реа) чаи)", 


я а). 1-8 


который полностью соответствует (29) при 5 =п- 7. 
Пример 5. В условиях примера 1 имеют место следующие соотношения: 


ЕР, =а-Ь К, =К,,=—а "а'. (30) 


Отсюда 


Согласно (26) 
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используя обратное преобразование, получим окончательный результат 
ы а-| _ 
ка, 
а 


который полностью соответствует (30) при $ =и-7. 


Для оценки полученных результатов по обобщению решений указанных типов 
дискретных аналогов интегральных уравнений Вольтерра П рода и исследованию 
новых их решений использованы программные средства пакета Зутбойс Май Тооох 
Майаь 6.5, который базируется на применении ядра символьной математики системы 
компьютерной алгебры Мар У В4 [9]. Приведенные выше утверждения и примеры про- 
верены с использованием языка символической математики и полностью подтверждены. 


Заключение 


Предлагаемый подход к установлению решений дискретных аналогов некоторых 
типов интегральных уравнений Вольтерра П рода на основе отыскания решений соот- 
ветствующих уравнений, связывающих ядро и резольвенту, дает возможность иссле- 
довать новые классы решений таких уравнений и упростить алгоритмы численной 
реализации при таблично заданных входных данных, что определяет существенные 
преимущества предлагаемых моделей при решении прикладных задач. Перспективы 
дальнейших исследований связаны с расширением круга возможных типов интеграль- 
ных уравнений Вольтерра П рода и их дискретных аналогов, для которых могут быть 
получены аналитические решения. 
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В.Ф. Миргород, 1.М. Гвоздева 

Дискретн! аналоги 1нтегральних рвнянь Вольтера П роду 

У статт! пропонуеться клас математичних моделей процесв у лиййних системах з факсованими станами. 
Так! математичн! модел е дискретними аналогами 1нтегральних равнянь Вольтера П роду. Розглянут! 
методи аналитичного розв’язання запропонованих равнянь на основ! вдшукання ришень вдповдних 
равнянь, що пов’язують ядро та резольвенту. Розглянута та доведена теорема, що встановлюе анаттичний 
вигляд резольвенти за заданим ядром для низки важливих окремих випадкв, а саме при сепарабельному 
вигляд] ядра. Розглянут! екввалентн! перетворення пропонованих математичних моделей до лйних 
разницевих равнянь. 


Г.Е. Мизогоа, Г.М. Суоз4еуа 

'Тве О5сгее Апаюгие$ оЁ УоКегга Пцеота! Едиайоп$ оЁ е 5есопа Кша 

ТБе с1а$ оР тафетайса] плоде! оЁ ргосез5ез ш фе Ппеаг зузветл$ \/И Ше Яхеа зваиез 15 оНеге шт Фе рарег. 
Зисн таетайса! ппо4е]$ аге Ше 415сгее апаюгиез ог Уо[егга ищеэга| едиайопз$ оЁ Фе зесопа Ка. ТБе 
тефо4$ о# апауНса1 зо[аноп$ оЁ Фе оЙеге4 едааНоп$ оп Ше Баз1$ оЁ Нп4ше оЁ зо[аНоп$ оЁ соггезроп те 
едфлайот$ Буш те Фе Кегпе! ап4 гезо[уеп аге сопз14егеа. ТБе еогет, дебпте Фе апа[уйса| Ка о гезоуепе 
оп Фе э1туеп Кегпе! ог зоте ппропапе зресла| сазез, патеЙу ай фе зерагае {уре оЁ Кегпе|, 15 сопз14егеа 
ап4 ргоуе4. Едилуа]ет{ тапзРогтаНоп$ оЁ 1е оЁеге4 тафетайса| по4е[$ © Фе Ппеаг Чегепсе едиаНоп$ 
аге сопз1дегед. 
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